Rocuronium pharmacokinetics and pharmacodynamics in the adductor pollicis and masseter muscles.

BACKGROUND: The aim of this study was to characterize the dose-effect relationship of rocuronium at the adductor pollicis and masseter muscles. METHODS: Ten, ASA I, adult patients, received a bolus dose of rocuronium 0.3 mg/kg during propofol based anesthesia. Train-of-four (TOF) was simultaneously monitored at the masseter and the adductor pollicis muscles until recovery. Rocuronium arterial serum concentrations were measured during 120 min. The first twitch of the TOF response was used to characterize the time-effect profile of both muscles using pharmacokinetic-pharmacodynamic analysis in NONMEM. A decrease in NONMEM objective function (∆OFV) of 3.84 points for an added parameter was considered significant at the 0.05 level. RESULTS: Onset time at the masseter (mean ± SD, 1.5 ± 0.9 min) was faster than at the adductor pollicis (2.7 ± 1.4 min, P < 0.05). Recovery, measured as the time to TOF ratio = 0.9 was similar between muscles 29.9 ± 6.7 (adductor pollicis) vs. 29.3 ± 8.1 (masseter). (P = 0.77). The estimated pharmacodynamic parameters [mean (95% CI)] of the adductor pollicis muscle and the masseter muscle were; plasma effect-site equilibration half-time (teq) 3.25 (2.34, 3.69) min vs. 2.86 (1.83, 3.29) min, (∆OFV 383.665); Ce50 of 1.24 (1.13, 1.56) mg/l vs. 1.19 (1.00, 1.21) mg/l, (∆OFV 184.284); Hill coefficient of 3.97 (3.82, 5.62) vs. 4.68 (3.83, 5.71), (∆OFV 78.906). CONCLUSIONS: We found that the masseter muscle has faster onset of blockade and similar recovery profile than adductor pollicis muscle. These findings were best, explained by a faster plasma effect-site equilibration of the masseter muscle to rocuronium.

A review of mixed-effects models of tumor growth and effects of anticancer drug treatment used in population analysis.

Population modeling of tumor size dynamics has recently emerged as an important tool in pharmacometric research. A series of new mixed-effects models have been reported recently, and we present herein a synthetic view of models with published mathematical equations aimed at describing the dynamics of tumor size in cancer patients following anticancer drug treatment. This selection of models will constitute the basis for the Drug Disease Model Resources (DDMoRe) repository for models on oncology.

Pharmacometrics: opportunity for reducing disease burden in the developing world: the case of Africa.

Pharmacometricians are virtually nonexistent in Africa and the developing world. The unrelenting burden of infectious diseases, which are often treated using medicines with narrow effectiveness and safety dose ranges, and the growing prevalence and recognition of non-communicable diseases represent significant threats for the patients, although affording an opportunity for advancing science. This article outlines the case for pharmacometricians to redirect their expertise to focus on the disease burden affecting the developing world.CPT: Pharmacometrics & Systems Pharmacology (2013) 2, e69; doi:10.1038/psp.2013.45; published online 28 August 2013.

Influence of obesity on propofol pharmacokinetics: derivation of a pharmacokinetic model.

BACKGROUND: The objective of this study was to develop a pharmacokinetic (PK) model to characterize the influence of obesity on propofol PK parameters. METHODS: Nineteen obese ASA II patients undergoing bariatric surgery were studied. Patients received propofol 2 mg kg(-1) bolus dose followed by a 5-20-40-120 min, 10-8-6-5 mg kg(-1) h(-1) infusion. Arterial blood samples were withdrawn at 1, 3, 5 min after induction, every 10-20 min during propofol infusion, and every 10-30 min for 2 h after stopping the propofol infusion. Arterial samples were processed by high-performance liquid chromatography. Time-concentration data profiles from this study were pooled with data from two other propofol PK studies available at http://www.opentci.org. Population PK modelling was performed using non-linear mixed effects model. RESULTS: The study involved 19 obese adults who contributed 163 observations. The pooled analysis involved 51 patients (weight 93 sd 24 kg, range 44-160 kg; age 46 sd 16 yr, range 25-81 yr; BMI 33 sd 9 kg m(-2), range 16-52 kg m(-2)). A three-compartment model was used to investigate propofol PK. An allometric size model using total body weight (TBW) was superior to all other models investigated (linear TBW, free fat mass, lean body weight, normal fat mass) for all clearance parameters. Variability in V2 and Q2 was reduced by a function showing a decrease in both parameters with age. CONCLUSIONS: We have derived a population PK model using obese and non-obese data to characterize propofol PK over a wide range of body weights. An allometric model using TBW as the size descriptor of volumes and clearances was superior to other size descriptors to characterize propofol PK in obese patients.

Morphine pharmacokinetics and pharmacodynamics in preterm and term neonates: secondary results from the NEOPAIN trial.

BACKGROUND: Relationships between plasma morphine concentrations and neonatal responses to endotracheal tube (ETT) suctioning are unknown in preterm neonates. METHODS: Ventilated preterm neonates (n=898) from 16 centres were randomly assigned to placebo (n=449) or morphine (n=449). After an i.v. loading dose (100 microg kg(-1)), morphine infusions [23-26 weeks postmenstrual age (PMA) 10 microg kg(-1) h(-1); 27-29 weeks 20 microg kg(-1) h(-1); and 30-32 weeks 30 microg kg(-1) h(-1)] were established for a maximum of 14 days. Open-label morphine (20-100 microg kg(-1)) was given for pain or agitation. Morphine assay and neonatal response to ETT suctioning was measured at 20-28 and 70-76 h after starting the drug infusion and at 10-14 h after discontinuation of the study drug. The concentration-effect response was investigated using non-linear mixed effects models. RESULTS: A total of 5119 data points (1598 measured morphine concentrations and 3521 effect measures) were available from 875 neonates for analysis. Clearance was 50% that of the mature value at 54.2 weeks PMA (CLmat(50)) and increased from 2.05 litre h(-1) 70 kg(-1) at 24 weeks PMA to 6.04 litre h(-1) 70 kg(-1) at 32 weeks PMA. The volume of distribution in preterm neonates was 190 litre 70 kg(-1) (CV 51%) and did not change with age. There was no relationship between morphine concentrations (range 0-440 microg litre(-1)) and heart rate changes associated with ETT suctioning or with the Premature Infant Pain Profile. CONCLUSIONS: A sigmoid curve describing maturation of morphine clearance is moved to the right in preterm neonates and volume of distribution is increased compared with term neonates. Morphine does not alter the neonatal response to ETT suctioning.

Mechanism-based concepts of size and maturity in pharmacokinetics.

Growth and development can be investigated using readily observable demographic factors such as weight and age. Size is the primary covariate and can be referenced to a 70-kg person with allometry using a coefficient of 0.75 for clearance and 1 for volume. The use of these coefficients is supported by fractal geometric concepts and observations from diverse areas in biology. Fat free mass (FFM) might be expected to do better than total body weight when there are wide variations in fat affecting body composition. Clearance pathways develop in the fetus before birth. The use of postnatal age as a descriptor of maturation is unsatisfactory because birth may occur prematurely; therefore postmenstrual age is a superior predictor of elimination function. A sigmoid E(max) model (Hill equation) describes gradual maturation of clearance in early life leading to a mature adult clearance achieved at a later age.

Developmental pharmacokinetics of morphine and its metabolites in neonates, infants and young children.

BACKGROUND: Descriptions of the pharmacokinetics and metabolism of morphine and its metabolites in young children are scant. Previous studies have not differentiated the effects of size from those related to age during infancy. METHODS: Postoperative children 0-3 yr old were given an intravenous loading dose of morphine hydrochloride (100 micro g kg(-1) in 2 min) followed by either an intravenous morphine infusion of 10 micro g h(-1) kg(-1) (n=92) or 3-hourly intravenous morphine boluses of 30 micro g kg(-1) (n=92). Additional morphine (5 micro g kg(-1)) every 10 min was given if the visual analogue (VAS, 0-10) pain score was >/=4. Arterial blood (1.4 ml) was sampled within 5 min of the loading dose and at 6, 12 and 24 h for morphine, morphine-3-glucuronide (M3G) and morphine-6-glucuronide (M6G). The disposition of morphine and formation clearances of morphine base to its glucuronide metabolites and their elimination clearances were estimated using non-linear mixed effects models. RESULTS: The analysis used 1856 concentration observations from 184 subjects. Population parameter estimates and their variability (%) for a one-compartment, first-order elimination model were as follows: volume of distribution 136 (59.3) litres, formation clearance to M3G 64.3 (58.8) litres h(-1), formation clearance to M6G 3.63 (82.2) litres h(-1), morphine clearance by other routes 3.12 litres h(-1) per 70 kg, elimination clearance of M3G 17.4 (43.0) litres h(-1), elimination clearance of M6G 5.8 (73.8) litres h(-1). All parameters are standardized to a 70 kg person using allometric 3/4 power models and reflect fully mature adult values. The volume of distribution increased exponentially with a maturation half-life of 26 days from 83 litres per 70 kg at birth; formation clearance to M3G and M6G increased with a maturation half-life of 88.3 days from 10.8 and 0.61 litres h(-1) per 70 kg respectively at birth. Metabolite formation decreased with increased serum bilirubin concentration. Metabolite clearance increased with age (maturation half-life 129 days), and appeared to be similar to that described for glomerular filtration rate maturation in infants. CONCLUSION: M3G is the predominant metabolite of morphine in young children and total body morphine clearance is 80% that of adult values by 6 months. A mean steady-state serum concentration of 10 ng ml(-1) can be achieved in children after non-cardiac surgery in an intensive care unit with a morphine hydrochloride infusion of 5 micro g h(-1) kg(-1) at birth (term neonates), 8.5 micro g h(-1) kg(-1) at 1 month, 13.5 micro g h(-1) kg(-1) at 3 months and 18 micro g h(-1) kg(-1) at 1 year and 16 micro g h(-1) kg(-1) for 1- to 3-yr-old children.

Population PKPD modelling of the long-term hypoglycaemic effect of gliclazide given as a once-a-day modified release (MR) formulation.

AIMS: To study the relationship between the pharmacokinetics (PK) of gliclazide and its long-term pharmacodynamic (PD) effect in a large population of Type 2 diabetic patients and to identify factors predicting intersubject variability. METHODS: A PKPD database of 634 Type 2 diabetic patients with a total of 5,258 fasting plasma glucose (FPG) samples was built up from the data collected during the clinical development of a modified release formulation of gliclazide (gliclazide MR). The PKPD analysis used a nonlinear mixed effect modelling approach. A mixture model was used to identify patients with a FPG response to treatment. In patients identified as responders, the decrease in FPG was related to gliclazide exposure (AUC) by an Emax relationship. An effect compartment was used to describe the link between PK and PD. A linear disease-progression model was used to assess the glycaemic deterioration observable over several months of treatment. Simulations were performed to evaluate the predictive performance of the PKPD model and to illustrate the time course of the antidiabetic effect of gliclazide MR. RESULTS: Disease state was found to be the main explanatory factor for intersubject variability in response to gliclazide. The percentage of responders to gliclazide, used as monotherapy, increased inversely to the number of classes of antidiabetic agents received prior to entry in the studies. In responders, the initial dose (30 mg) of the gliclazide MR dosing regimen induced half of the maximum hypoglycaemic effect. The equilibration half-life between the PK and PD steady states was 3 weeks (intersubject variability of 84%). The rate of disease progression was 0.84 mmol l(-1) year(-1) (intersubject variability 143%). The PKPD model adequately predicted the FPG profiles of 234 patients who received the current formulation of gliclazide. Simulation of a 1-year parallel dose ranging clinical trial illustrated the influence of dose, time and type of previous antidiabetic treatment on the percentage of patients with clinically significant improvement of blood glucose control. CONCLUSIONS: This population PKPD analysis has characterized the relationship between the exposure to gliclazide and its long-term hypoglycaemic effect, and has established that the intersubject variability in response is mostly related to disease state. These results underline the clinical interest of quickly increasing the dose of gliclazide MR according to the response to treatment in order to achieve effective blood glucose control.

A pharmacokinetic model to predict the PK interaction of L-dopa and benserazide in rats.

PURPOSE: To study the PK interaction of L-dopa/benserazide in rats. METHODS: Male rats received a single oral dose of 80 mg/kg L-dopa or 20 mg/kg benserazide or 80/20 mg/kg L-dopa/benserazide. Based on plasma concentrations the kinetics of L-dopa, 3-O-methyldopa (3-OMD), benserazide, and its metabolite Ro 04-5127 were characterized by noncompartmental analysis and a compartmental model where total L-dopa clearance was the sum of the clearances mediated by amino-acid-decarboxylase (AADC), catechol-O-methyltransferase and other enzymes. In the model Ro 04-5127 inhibited competitively the L-dopa clearance by AADC. RESULTS: The coadministration of L-dopa/benserazide resulted in a major increase in systemic exposure to L-dopa and 3-OMD and a decrease in L-dopa clearance. The compartmental model allowed an adequate description of the observed L-dopa and 3-OMD concentrations in the absence and presence of benserazide. It had an advantage over noncompartmental analysis because it could describe the temporal change of inhibition and recovery of AADC. CONCLUSIONS: Our study is the first investigation where the kinetics of benserazide and Ro 04-5127 have been described by a compartmental model. The L-dopa/benserazide model allowed a mechanism-based view of the L-dopa/benserazide interaction and supports the hypothesis that Ro 04-5127 is the primary active metabolite of benserazide.

The pharmacokinetics of remifentanil in patients undergoing coronary artery bypass grafting with cardiopulmonary bypass.

UNLABELLED: Remifentanil is a potent opioid with a short duration of action. It has the potential for large-dose opioid anesthesia without an obligatory prolonged period of mechanical ventilation. However, because of high clearance and rapid tissue distribution, cardiopulmonary bypass (CPB) may influence its pharmacokinetics and alter drug requirements. We administered remifentanil by continuous infusion to 68 patients having coronary artery bypass graft surgery during CPB with hypothermia to describe the effects of these interventions on its pharmacokinetics. Remifentanil concentrations were measured before, during, and after CPB. Disposition was best described by a two-compartment model. The volume of distribution increased by 86% with institution of CPB and remained increased after CPB. Elimination clearance decreased by 6.37% for each degree Celsius decrease from 37 degrees C. IMPLICATIONS: Remifentanil concentrations decrease with the institution of cardiopulmonary bypass because of an increase in the volume of distribution. The decrease in elimination clearance with hypothermia results in increased total remifentanil concentrations during cardiopulmonary bypass if the infusion rate is not altered. More constant blood remifentanil levels may be obtained by reducing remifentanil infusion rate by 30% for each 5 degrees C decrease in temperature.

Target concentration intervention: beyond Y2K.

Target concentration intervention (TCI) is proposed as an alternative conceptual strategy to therapeutic drug monitoring (TDM). It is argued that the idea of a therapeutic range has limited the interpretation of measured drug concentrations and diminished the anticipated clinical benefit to patients by use of an oversimplified pharmacodynamic model. TCI on the other hand embraces pharmacokinetic and pharmacodynamic concepts and uses the idea of a target effect and associated target concentration to make rational individual dose decisions.

Drug treatment effects on disease progression.

Degenerative diseases are characterized by a worsening of disease status over time. The rate of deterioration is determined by the natural rate of progression of the disease and by the effect of drug treatments. A goal of drug treatment is to slow disease progression. Drug treatments can be categorized as symptomatic or protective. Symptomatic treatments do not affect the rate of disease progression whereas protective treatments have the ability to slow disease progression down. Many current methods for describing disease progression have two common drawbacks: a linear relationship between time and disease status is assumed, and within- and between-subject variability is ignored. Disease progress models combined with pharmacokinetic pharmacodynamic models and hierarchical random effects statistical models provide insights into understanding the time course and management of degenerative disease.

Acetaminophen analgesia in children: placebo effect and pain resolution after tonsillectomy.

BACKGROUND: Pharmacodynamic models of acetaminophen analgesia in children have not explored the efficacy of single oral doses greater than 40 mg/kg. METHODS: Children aged 9.0 +/- 3.0 years (+/- SD) and weight 37.9+/- 16.6 kg undergoing outpatient tonsillectomy were randomised to receive acetaminophen elixir 40 mg/kg (n = 12). high dose acetaminophen elixir 100 mg/kg (n =20) or placebo (n=30) 0.5 -1 h preoperatively. No other analgesics were given. Individual acetaminophen serum concentrations and pain scores [visual analogue scale (VAS) 0-10] were measured over a 4-8 h postoperative period. These data were pooled with data from a previous study investigating acetaminophen pharmacodynamics (n = 120) and analysed using a non-linear mixed effect model. Placebo effects and drug effects were modelled using effect-site concentration models. RESULTS: A one-compartment model with first-order input, lag time and first-order elimination was used to describe the population pharmacokinetics of acetaminophen. Pharmacokinetic parameter estimates were similar to those previously described. Pharmacodynamic population parameter estimates [population variability coefficient of variation (CV)] for a maximum analgesic effect (Emax) model, in which the greatest possible pain relief (VAS 0-10) equates to an Emax of 10, were Emax 5.17 (64%) and 50% effective concentration 9.98 mg/l (107%). The equilibration half-life (t(eq)) of the analgesic effect compartment was 53 min (217%). A placebo drug model for the effects of placebo response had a t(eq) of 1.96 h (40%), an elimination half-life of 2.06 h (50%) and a potency of 1.54 pain relief units (24%). CONCLUSIONS: High dose acetaminophen (100 mg/kg) was no more effective than 40 mg/kg and was associated with increased nausea and vomiting. A target effect compartment concentration of 10 mg/l is expected to produce a pain reduction of 2.6 units. The placebo model accounted for a maximum pain reduction of 5.6 units at 3 h. The combination of placebo effect and preoperative acetaminophen 40 mg/kg results in pain scores below 4 units for 5 h postoperatively.

Prediction of the outcome of a phase 3 clinical trial of an antischizophrenic agent (quetiapine fumarate) by simulation with a population pharmacokinetic and pharmacodynamic model.

A completed phase 3 trial result was simulated 100 times on the basis of a simulation model of quetiapine fumarate (Seroquel), an antischizophrenic agent. The simulation was executed by analysts who were completely blinded from results of the actual trial until after the simulations were submitted to the holder of the trial results. Data from two clinical investigations of quetiapine in patients with schizophrenia were analyzed by use of nonlinear mixed effects modeling to derive a population pharmacokinetic- and pharmacodynamic-based simulation model. The time course of quetiapine concentrations was described by use of a one-compartment open linear pharmacokinetic model with first-order absorption and elimination. The combination of an inhibitory maximum effect pharmacodynamic model for the active treatment effect and a linear function of time for the placebo effect characterized the observed time course of change in the Brief Psychiatric Rating Scale. Simulation results were compared with those in the actual trial to evaluate how well the simulations predicted the outcome. The actual trial results for all doses except the placebo group fell within the predicted Brief Psychiatric Rating Scale scores +/- 1 SE. Unlike the phase 2 trial, from which the pharmacokinetic/pharmacodynamic model was developed, the placebo group in the actual phase 3 trial showed deterioration of Brief Psychiatric Rating Scale scores with time. We conclude that variable placebo responses observed in short-term studies of schizophrenia provide an inadequate basis for the modeling and simulation of placebo subjects in clinical trials. Knowledge of the range of placebo response observed in other studies may have provided an improved basis for the placebo effect model. The model for active drug produced adequate predictions of the actual trial outcomes.

Investigations using logistic regression models on the effect of the LMA on morphine induced vomiting after tonsillectomy.

The effect of intraoperative airway management on postoperative vomiting after tonsillectomy is unknown. Logistic regression was used in a retrospective study to investigate the effect of the laryngeal mask airway (LMA) on a morphine dose-vomiting response curve. Charts were reviewed in 351 children in whom the airway was managed with either a LMA (n=177) or a tracheal tube (n=174). A mean perioperative morphine dose of 0.10 mg.kg(-1) (SD 0.09) was used in 248 children and a further 103 children were given no opioid. One hundred and eighteen of these 248 children vomited (47.6%) compared to 14 of 103 children given no morphine (13.6%). The probability of vomiting was related to morphine dose using logistic regression with both a linear and an E(max) model. Both the calibration (Hosmer-Lemishow goodness of fit chi-squared test lambda(2), P=0.81) and discrimination (area under the receiver operating characteristic plot, AUC ROC=0.67) of the E(max) model were better than the linear model (lambda(2), P=0.49; AUC ROC=0.64). Pharmacodynamic parameter estimates for the Emax model were P(0) (the baseline probability of vomiting) 0.139, P(max) (the maximal probability of vomiting due to morphine) 0.96, ED(50) (morphine dose that induces an effect equivalent to 50% of the logit P(max)) 0.09 mg.kg-1. The probability of vomiting was 50% after morphine 0.125 mg.kg-1. The use of the LMA had no effect on this dose-response curve. A covariate analysis investigating propofol for induction or isoflurane for the intraoperative maintenance of anaesthesia, however, showed that both drugs shifted the curve to the right. The probability of vomiting was 50% after morphine 0.17 mg.kg(-1) and 0.21 mg.kg(-1) for the isoflurane and propofol use curves, respectively. The concomitant use of propofol and isoflurane, but not the use of the LMA, decreases the probability of vomiting due to morphine.

A model for size and age changes in the pharmacokinetics of paracetamol in neonates, infants and children.

AIMS: The aims of this study were to describe paracetamol pharmacokinetics in neonates and infants. METHODS: Infants in their first 3 months of life (n = 30) were randomised to sequentially receive one of three paracetamol formulations (dose 30-40 mg kg-1) over a 2 day period. The formulations were (a) elixir, (b) glycogelatin capsule suppository and (c) triglyceride base suppository. Approximately six blood samples were taken after each dose over the subsequent 10-16 h. Data were analysed using a nonlinear mixed effect model. These neonatal and infant data were then included with data from four published studies of paracetamol pharmacokinetics (n = 221) and age-related pharmacokinetic changes investigated. RESULTS: Population pharmacokinetic parameter estimates and their coefficients of variation (CV%) for a one compartment model with first order input, lag time and first order elimination were volume of distribution 69.9 (18%) l and clearance 13.0 (41%) l h-1 (standardized to a 70 kg person). The volume of distribution decreased exponentially with a half-life of 1.9 days from 120 l 70 kg-1 at birth to 69.9 l 70 kg-1 by 14 days. Clearance increased from birth (4.9 l h-1 70 kg-1) with a half-life of 3.25 months to reach 12.4 l h-1 70 kg-1 by 12 months. The absorption half-life (tabs) for the oral preparation was 0.13 (154%) h with a lag time (tlag) of 0.39 h (31%). Absorption parameters for the triglyceride base and capsule suppositories were tabs 1.34 (90%) h, tlag 0.14 h (31%) and tabs 0.65 (63%) h, tlag 0.54 h (31%), respectively. The tabs for elixir and capsule suppository in children under 3 months were 3.68 and 1.51 times greater than children over 3 months. The relative bioavailability of rectal formulations compared with elixir were 0.67 (30%) and 0.61 (23%) for the triglyceride base and capsule suppositories, respectively. CONCLUSIONS: Total body clearance of paracetamol at birth is 62% and volume of distribution 174% that of older children. A target concentration above 10 mg l-1 in approximately 50% subjects can be achieved by a dose from 45 mg kg-1 day-1 at birth and up to 90 mg kg-1 day-1 in 5-year-old children. A reduced dose of 75 mg kg-1 day-1 in an 8-year-old child is sufficient because clearance is a nonlinear function of weight.

Simulation of clinical trials.

Computer simulation of clinical trials has evolved over the past two decades from a simple instructive game to "full" simulation models yielding pharmacologically sound, realistic trial outcomes. The need to make drug development more efficient and informative and the awareness that many industries make extensive use of simulation in product development have advanced considerably the use of simulation of clinical trials in pharmaceutical product development over the past decade. The structural and stochastic components of trial simulation models are explained as a prelude to a listing of representative simulation projects, reflecting investigative applications of statistical methods, trial design comparisons, and full simulation of new drugs being developed. Lessons learned from these projects are reviewed in the context of their current impact and potential for influencing the future of drug development.

The dose-effect relationship for morphine and vomiting after day-stay tonsillectomy in children.

A dose-response curve for intravenous morphine and vomiting was investigated in children having day-stay tonsillectomy. A retrospective chart review was performed for the 164 children fulfilling the inclusion criteria. Morphine (mean 0.09 mg/kg SD 0.05) was used in 108 children in the perioperative period and a further 56 children were given no opioid. Fifty-five of these 164 children vomited and 20 children required an overnight stay in hospital because of vomiting. The probability of vomiting or overnight stay in hospital was related to morphine dose (by logistic regression). The overall probability of vomiting after morphine 0.1 mg/kg was 50% and the probability of admission for vomiting with this dose was 10%. Pharmacodynamic parameter estimates for postoperative vomiting were P0 (the baseline probability of vomiting, with no opioid) 0.115, Pmax (the maximal probability of vomiting due to morphine) 0.997, ED50 (morphine dose that induces an effect equivalent to 50% of the logit Pmax) 0.18 mg/kg. Parameter estimates for overnight stay because of vomiting after morphine administration were P0 0.038, Pmax 0.999, ED50 0.369 mg/kg. Satisfactory postoperative analgesia in children has been reported with morphine 0.05 to 0.15 mg/kg. Doses above 0.1 mg/kg are associated with a greater than 50% incidence of vomiting. Our data suggests that lower doses of morphine are associated with a decreased incidence of emesis after tonsillectomy in children.